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Introduction 

This  report  describes  a  computational  routine  to 
figure  average  skid  distance  for  an  irregular  area  with 
or  without  variable  log  density  and  presents  a  method 
to  compute  average  yarding2  distance  with  a  hand- 
held calculator.  Theory  underlying  this  computer 
program  is  presented  in  an  appendix,  as  is  theory 
underlying  the  application  of  the  variable  log  density 
formula,  but  understanding  of  this  theory  is  not 
necessary  to  use  the  computer  program. 

Methods  currently  available  for  computing  average 
skid  distance  assume  that  logs  are  evenly  distributed 
over  the  area.  While  this  assumption  may  have  been 
satisfactory  in  the  past,  current  timber  harvest  prac- 
tice can  cause  different  log  densities  on  subareas  of 
the  same  skid  area.  Average  skid  distance  can  be 
significantly  different  if  variable  log  density  is  con- 
sidered by  subarea. 

Past  Work 

Average  skid  distance  has  been  computed  for 
many  years  in  sale  appraisal  and  logging  planning. 
Matthews  (1942)  used  the  method  of  equal  areas  to 
derive  average  skid  distance  (ASD)  factors  for  areas 
shaped  like  circles,  wedges,  squares,  and  rectangles. 
Other  investigators  have  shown  that  Matthews'  fac- 
tors for  the  circle  and  wedge  are  in  error  (Lysons  and 
Mann  1965).  They  show  that  average  skid  distance 
for  these  two  shapes  is  two-thirds  of  the  external 
(maximum  possible)  skidding  distance  (ESD),  not 
0.707  of  the  external  skidding  distance  as  Matthews 
computed. 

Suddarth  and  Herrick  (1964)  computed  formulas 
based  on  integral  calculus  for  average  skid  distance 
of  areas  shaped  like  circles,  wedges,  right  triangles, 


2  In  this  paper,  the  words  yarding  and  skidding  are  used 
interchangeably.  The  computing  methods  presented  apply 
either  to  cable  yarding  or  tractor  skidding  distance  computa- 
tion. Distance  computation  for  cable  yarding,  however,  does 
not  include  lateral  yarding  distance  from  the  cable  road. 


and  rectangles.  Based  on  their  work,  Matthews'  skid 
distance  factor  for  a  rectangle  is  also  questionable. 
According  to  Matthews,  ASD  =  0.578  (ESD),  but 
Suddarth's  formula  gives  the  result,  ASD  =  0.593 
(ESD).  This  is  only  a  2  Vi%  error,  which  is  negligible 
for  most  logging  applications  since  other  factors 
affecting  skid  production  introduce  more  serious 
variations.  However,  the  difference  should  be  noted. 

Unfortunately,  skidding  areas  are  not  always  shaped 
like  circle  segments  or  rectangles.  Integral  calculus 
methods  using  numerical  approximations  for  com- 
puting average  skid  distance  have  been  extended  to 
irregularly  shaped  areas  (Peters  and  Burke  1972).  A 
desk  top  computer  and  digitizer  facilitate  computa- 
tion of  average  skid  distance  and  total  area  within  the 
prospective  skidding  area.  However,  a  digitizer  is  not 
always  readily  available,  so,  the  numerical  method 
described  here  is  an  efficient  substitute  which  can  be 
applied  on  a  hand-held  programable  calculator.  If 
this  equipment  also  is  unavailable,  computations  can 
be  done  on  a  desk-top  calculator  with  storage  reg- 
isters, although  analysis  time  will  then  increase. 
Appendix  1  presents  basic  program  logic  and  infor- 
mation. 


Different  Log  Densities 

A  yarding  area  may  have  two  or  more  distinctly 
different  log  densities.  How  this  condition  affects 
average  yarding  distance  (AYD)  depends  on  the 
density  difference  and  the  location  of  each  different 
log  density  area  with  respect  to  the  landing. 

A  region  of  heavy  log  density  located  far  from  the 
landing  will  significantly  increase  average  yarding 
distance.  The  opposite  is  also  true.  How  much 
difference  in  skid  distance  is  involved  in  these  two 
extremes  can  be  determined  using  the  methods 
described  here. 

There  are  actually  two  distinct  parts  to  the  method 
presented.  The  first,  the  computer  program,  helps 
compute  average  yarding  distance  for  a  sub-area 
where  log  density  is  assumed  to  be  uniform.  The 
method  is  not  radically  different  from  those  cited 
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earlier,  except  for  its  theoretical  base,  and  no  digiti- 
zer is  required.  The  second  part  allows  the  logging 
analyst  to  combine  two  or  more  areas  of  different  log 
densities  to  find  an  overall  average  yarding  distance. 


Irregular  Boundaries  and  Different  Log  Densities 

First,  log  density  is  assumed  to  be  uniform  over  the 
entire  area.  Figures  for  total  area  and  average  skid 
distance  of  the  skidding  area  shown  in  figure  1  are 
computed  by  the  hand-held  caluclator  and  compared 
to  results  from  the  digitizer-computer  program: 


Total  area 

Aver,  skid  distance 


Computer- 
digitizer 

50.00  acre 
1,012  feet 


Hand-held 
calculator 

50.51  acre 
1,021  feet 


Next,  consider  the  three  areas  with  different  log 
density  shown  in  figure  1,  and  compute  acreage  and 


average  skid  distance  (ASD)  for  each  area  with  the 
programable  calculator: 


Area,  i 

1 

2 

3 


Area  log 
density,  8. 
logs/ acre 

400 
100 
200 


Area,  A.     ASD,  r. 


acres 

17.25 
22.87 
10.39 


feet 

1,466.0 
697.7 
992.3 


Average  skid  distance  for  the  entire  cutting  unit 
with  three  areas  of  differing  log  density  is  1,222.6  feet 
(appendix  2). 

For  the  entire  area,  neglecting  different  log  densi- 
ties, computed  results  from  the  digitizer-computer 
method  and  the  hand-held  programable  calculator 
are  similar.  When  the  area  is  divided  into  three  log 
density  areas  with  the  highest  log  density  far  from  the 
landing  in  area  1,  the  average  skid  distance  increases 
about  200  feet,  showing  the  effect  of  considering 
different  log  densities. 


-500-1- 


Figure  1. -Irregular  boundary  skidding  area  with  three  subareas  of  differing  log  density. 
Coordinate  values  are  in  feet.  (Based  on  Peters  and  Burke  1972,  p.  8) 
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In  contrast,  suppose  the  log  density  in  area  2 
nearest  the  landing  is  400  logs/acre  and  the  density  in 
area  1  far  from  the  landing  is  100  logs/acre,  the 
average  skid  distance  would  then  be  847  feet,  about 
160  feet  shorter  than  average  skid  distance  computed 
without  regard  for  log  density. 

No  effect  due  to  slope  was  considered  in  the 
previous  computations.  If  the  skidding  area  is  gener- 
ally level,  these  skid  distances  would  be  approxi- 
mately correct.  However,  when  average  overall  slope 
is  32%  or  greater,  the  actual  skid  distance  on  the 
ground  is  5%  or  larger  than  the  computed  horizontal 
distance  neglecting  slope.  These  figures  give  some 
indication  of  slope  effects  on  accuracy.  Since  most 
tractor  logging  is  done  on  slopes  less  than  32%, 
computed  average  skid  distance  will  be  accurate 
within  5%  of  the  value  for  flat  ground  in  these 
instances.  For  any  slope,  but  especially  those  greater 
than  32%,  average  yarding  distance  can  be  corrected 
for  slope  as  follows: 


ASD  (corrected)  = 


ASP  (horizontal) 


cos  [tan  -'(percent  slope/100] 


Another  use  for  the  methods  described  here  is  to 
compute  the  components  of  average  skid  distance 
where  logs  are  skidded  to  one  or  more  in-woods  swing 
decks  before  being  forwarded  to  the  main  concentra- 
tion landing.  Average  skid  distance  to  each  respec- 
tive swing  deck  is  computed  just  as  previously 
described.  If  logs  are  uniformly  distributed  only  area 
is  considered;  but  if  one  or  more  separate  log 
densities  are  found,  then  both  area  and  variable  log 
density  are  taken  into  account.  Average  skid  distance 
from  a  swing  deck  to  the  central  landing  is  the  same  as 
the  distance  between  the  two  points  because  all  logs 
are  concentrated  at  the  swing  landing  before  for- 
warding. 


Discussion 

Logs  Carried  Per  Turn 

One  of  the  assumptions  underlying  the  previous 
analysis  is  that  the  average  number  of  logs  carried  on 
each  turn  can  be  defined.  In  the  derivation  in 
appendix  2,  this  assumption  allows  the  expression  for 
average  skid  distance  to  be  independent  of  number  of 
logs  carried  per  turn.  This  assumption  is  somewhat 
easier  to  justify  for  cable  choker  skidders  than  for 
grapple  skidders. 


For  cable  skidders,  average  logs  per  turn  is  a 
function  of  number  of  chokers  carried,  of  average  log 
size,  of  log  dispersion  and  orientation  with  respect  to 
the  skid  trail,  and  of  operator  choice.  For  a  given  area 
where  these  variables  can  be  defined,  average  logs 
per  turn  can  also  be  estimated. 

For  grapple  skidders,  all  factors  above,  except 
number  of  chokers,  determine  number  of  logs  carried 
per  turn.  In  addition,  grapple  size  will  affect  the 
number  of  logs  carried. 

The  analysis  in  appendix  2  is  also  valid  if  weight 
per  acre,  wai  or  volume  per  acre,  vgi,  is  substituted  for 
logs  per  acre,  For  tractor  skidding,  especially 
grapple  skidders,  and  for  cable  yarding,  weight  is  the 
limiting  factor  on  load  size  for  a  given  set  of  circum- 
stances. Weight  can  be  approached  indirectly  by 
using  the  species  weight  per  cubic  foot  and  cubic 
volume  per  acre  in  the  derivations  of  appendix  2,  and 
the  computer  procedures  of  appendix  1. 

Wander  Factor 

Tracked  or  rubber  tired  skidders  seldom  travel  in  a 
straight  line  from  load  point  to  deck.  Actual  distance 
traveled  is  usually  greater  than  the  computed  straight- 
line  average  skid  distance.  Thus,  the  computed 
average  skid  distance  can  be  multiplied  by  a  wander 
factor  3  to  better  estimate  distance  skidded.  However, 
the  wander  factor  will  not  change  the  effect  of 
variable  log  density.  If  an  area  of  high  log  density  is 
located  far  from  the  deck,  average  skid  distance  will 
be  greater  than  if  logs  were  uniformly  distributed 
over  the  entire  skid  area.  If  a  wander  factor  is  not 
known,  then  the  best  that  can  be  done  is  to  compute  a 
straight-line  skid  distance,  knowing  that  actual  skid 
distance  will  be  greater  by  some  unknown  amount. 


Application 

The  average  skid  distance  computed  by  the  method 
described  here  is  simply  the  average  one-way  distance, 
weighted  by  number  of  logs  per  area,  that  is  traveled 
by  a  skidder  or  cable  yarder  carriage.  Not  included  in 
the  computation  are  distances  traveled  to  attach 
chokers  from  a  skidder  equipped  with  a  winch  cable 
or  lateral  distances  traveled  to  attach  chokers  from  a 
cable  yarder  equipped  with  a  skid  line. 

A  survey  of  timber  sale  appraisal  techniques  in- 
dicates that  external  and  average  yarding  distance 


3  Wander  factor  is  the  ratio  of  actual  travel  distance  to 
straight  line  distance. 
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are  key  factors  in  computing  estimated  yarding  cost. 
Lateral  yarding  distance,  however,  is  not  generally 
included  with  average  yarding  distance,  and  con- 
sequently, was  not  used  in  this  analysis. 

Summary 

When  individual  landings  and  the  area  each  serves 
can  be  identified  as  part  of  timber  sale  planning, 
methods  presented  here  can  help  compute  average 
skid  distance.  In  addition,  areas  of  different  known 
log  densities  can  be  included  in  the  analysis.  Effects 
of  differing  log  densities  on  average  skid  distance  can 
be  significant. 

The  interested  user  can  combine  these  examples 
with  theory  presented  in  the  appendices  to  alter  the 
method  for  his  own  application.  These  computations 
may  also  be  programed  on  other  computers,  including 
large  scale  digital  computers. 
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Appendix  1 


Average  Skid  Distance  Computations 

A  Numerical  Approach 

The  essential  feature  of  this  numerical  method  is 
that  an  area  with  any  shape  can  be  bordered  and 
approximated  by  a  series  of  connected  straight  line 
segments.  Lines  drawn  from  the  ends  of  each  such 


line  segment  to  the  coordinate  origin  form  a  triangle 
which  is  the  basic  area  for  subsequent  computation 
(fig.  2). 

Consider  for  the  moment,  only  the  single  triangle 
highlighted  in  figure  2.  Average  skid  distance  for  this 
triangle  is  the  distance,  ^  from  the  origin  of  coordi- 
nates, 0,  to  the  triangle's  polar  centroid,  P.C.,  assum- 
ing uniform  density  within  the  triangle. 


2000  h 


0^600,1800) 


Line  segment  3 
(x,  ,yj) 
(1600,1400) 

Mi, 

Line  segment  2  =  ^3 

(x,.,,yM  ) 

(1800,800) 

PC, Polar  Centroid 
of  the  triangle 


1 — r 
I  1 

Line  |secjn|ient 


(900,300) 
r,  Average  skid  dijatqnde  for  the  triangle 
 i      '  i  i  !  i  i  


1000 
Feet 


Xl  +  Xl-l 

~ 2- 


2000 


Figure  2.— Irregular  area  boundary  approxi- 
mated by  a  polygon.  The  shaded  triangle  is 
the  basic  area  for  computation  of  average 
skid  distance  (ASD)  and  area. 
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The  exact  formula  for  distance,  is 


f1   T  C2 


+ 


ioge  12  3 


/  +t  -e 

1  2  3 


entered  in  clockwise  sequence,  the  resulting  area  is 
positive.  The  absolute  value  function  in  the  calculator 
program  is  applied  to  the  final  value  for  area  and 
insures  positive  area  consistent  with  expectation. 

The  average  skid  distance  for  the  entire  area 
shown  in  figure  2  is 

n 

X  Ai  *i 

ASD  =  !__  i  =  1,  2,      no.  of  boundary 

y  line  segments. 


where 


'2  =  v2 


+y2, 
1   » 1-1 


'I  =x2  +  y2 

2  I        '  I 


and  ^  =  (xi-xi-i>2+(yi-yi-i)2 

are  the  lengths  of  the  triangle  sides,  and  where  the 
coordinates,  (xj,  and  (\v  y^),  represent  the  two 
non-zero  verticles  of  each  computing  triangle  as 
shown  in  figure  24.  Details  of  the  derivation  can  be 
obtained  from  the  author. 

An  alternative  derivation  is  based  on  approximating 
the  distance,     by  the  expression 


AYD=^r  ^2/3 


AYD=*^V3 


x:  +x 


i  -1 


+ 


V;  +y 


I  -1 


1/2 


(  xi  +xi-i  )2  +  )2 


V2 


This  approximate  formula  requires  one  vertex  of  the 
triangle  be  located  at  the  origin,  o.  The  approxi- 
mation improves  as  the  angle  of  the  triangle  touching 
the  origin  decreases. 

The  triangle  area,  Aj,  is  defined  by  coordinates  of 
the  line  segment  endpoints  as 

Aj  =V2(xiyi_1-xi_1  yj) 

The  sign  of  A.  and  of  the  sum  of  all  Aj's  depends  on 
the  sequence  of  data  entry.  For  example,  counter- 
clockwise data  entry  in  figure  2  results  in  an  area  with 
a  negative  sign.  This  is  not  an  error,  but  a  result  of 
data  entry  sequence.  If  data  for  this  example  were 

4  Personal  communication  from  P. A.  Peters,  Associate  Pro- 
fessor, Forest  Engineering  Department,  Oregon  State  Univ., 
Corvallis,  Oreg.,  1976. 


Table  1  illustrates  the  mechanics  of  computation  for 
average  skid  distance. 

As  stated  earlier,  the  minus  sign  in  the  column  total 
for  triangle  area  and  in  the  column  for  polar  moment 
results  from  the  counter-clockwise  direction  of  line 
segment  numbering  in  figure  2. 


Table  1. — Illustration  of  mechanical  computations  for  obtaining 
average  skid  distance 


Triangle 


Polar 


(line  segment)  coordinate 


Triangle 
area 


First  polar 
moment, 


number,  i 

r. 
i 

A 

r.A. 

1 

971.83 

-90 

X 

103 

-0.8746 

X 

108 

2 

1,349.90 

-620 

X 

103 

-8.3694 

X 

108 

3 

1,294.43 

-1,020 

X 

103 

-13.2032 

X 

108 

4 

954.52 

+300 

X 

103 

+2.8636 

X 

108 

5 

620.04 

+270 

X 

103 

+  1.6741 

X 

108 

Total 

-1,160 

X 

103 

-17.9095 

X 

108 

ASD=  -17.9095  X108 
-1,160x103 


Area-       -1,160  X  103  Sq.ft.  I 

Area  -  I  — — — —  —~  1  =  26.63  acres 

\\   43,560 sq.ft. /acre  I 


Program  for  a  Hand-Held  Computer 

The  preceding  algorithm  is  programed  on  the 
Hewlet-Packard  65,  a  hand-held  programable  cal- 
culator5. Figures  3  and  4  show  the  coding  and  opera- 
ting steps. 


5  The  use  of  trade,  firm,  or  corporation  names  in  this  publica- 
lication  is  not  an  official  endorsement  or  approval  by  the  U.S. 
Department  of  Agriculture  of  any  product  to  the  exclusion  of 
others  which  may  be  suitable. 
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HP-65  Program  Form 

Skid  distance  —  irregular  area 
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09 

+ 
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Rq 

f 

31 

l|(xi  +  xi.1)  +  (y;  +  yj.,)2]  14 

V  x 

09 

 O  

3 

03 

\  'A  [(xj  +  Xj.,)2  +  (y;  +  yj_,)*|  '/2 

LABELS 

^  Initialize 

 m — 

o  1 

30      STO  7 

33  07 

Store  r:  in  R7 

80 

g         Aver  skid  disL 

RCL  3 

34  03 

Recall  xj 

Q      Skid  Area 

RCL  2 

34  02 

Recall  yj_, 

D 

X 

71 

X,  >'i-l 

F 

RCL  1 

34  01 

Recall  Xj., 

O 

RCL  4 

34  04 

Recall  yj 

■J  Compute 

X 

71 

xi-i  y; 

2 

— 

51 

(xj  yj.,  -  Xj.,  yj) 

3 

2 

02  J 

}  A;  =  <A  (x;  y:  ,  •  X:.,  yj) 

I        1               v    1  J  11         11  J  V 

4 

81 

5 

40      STO  8 

33  08 

Store  Aj  in  R8 

90 

6 

RCL  5 

34  05 

Recall  current  value  £Aj 

7 

+ 

61 

Add  new  A;  value 

8 

STO  5 

33  05 

Store  new  cur.  value  lAj 

9 

RCL  7 

34  07 

Recall  rj 

Figure  3.— Program  coding  for  the  average  skid  dis- 
tance program  computed  on  an  HP-65  programmable 
calculator. 
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FLAGS 
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50  RCL  :! 

34  03 
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wo 
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HP-65  User  Instructions 


Title. 


Skid  distance  —  irregular  area 


Programmer , 


D.  Donnelly 


I  ASP   |  ^d"T 


Page. 
Date . 


.of 


1=3    In   I  i 


t — r 


STEP 


INSTRUCTIONS 


INPUT 
DATA/UNITS 


KEYS 


OUTPUT 
DATA/UNITS 


Key  first  x  coord,  xq;  Key/ENTERI/ 


Key  first  y  coord.,  y0;  Key/A7 


When  program  stops  w/xQ  in  display. 


Key  next  x  coord.,  x;;  Kev/R/S/ 


R/S 


Key  next  y  coord.,  yj;  Kev/R/S/ 


R/S 


Program  stops  with  last  x  coord,  entered 


display; 


If  more  coordinates  to  input,  go  to  #4 


When  closure  ^  is  achieved  and  processed, 


to  7 


Kev/RTN7./B7 


RTN 


Aver,  skid 
dist.,  ft 


Key  EJ 


Skid  area,  ac 


To  start  a  new  area,  Key  /17/STK/  ,  /f7/REG/ 


STK 


/RTN/ 


REG 


RTN 
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Go  to  Step  1 


Closure  is  achieved  when  all  coordinates  have  been  entered,  including,  at  the  end, 
the  values  ofxQ  &  yQ;  ie.,  xg  &  yg  are  entered  twice,  once  at  the  start  and  once  at  the  end. 


Figure  4.— Operating  instructions  for  the  average  skid 
distance  program  computed  on  an  HP-65  program- 


mable calculator. 
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Two  operational  points  should  be  noted.  The 
actual  periphery  of  an  area  is  more  closely  estimated 
by  increasing  the  number  of  approximating  line 
segments.  At  some  point,  though,  increased  accuracy 
from  additional  line  segments  is  offset  by  increased 
computational  time.  Second,  since  we  are  simulating 
a  continuous  curve  by  a  series  of  straight  line 
segments,  careful  placement  of  a  given  number  of 


line  segments  in  relation  to  the  area  boundary  can 
increase  accuracy.  Sharp  changes  in  the  boundary 
direction  are  good  places  to  start  or  end  a  straight 
segment. 

The  interested  reader  can  use  figure  1  as  an  exer- 
cise and  compute  the  results  previously  given  in 
Irregular  Boundaries  and  Different  Log  Densities 
using  the  program  code  and  operating  instructions. 


Appendix  2 
Variable  Log  Density 


One  Possible  Derivation 

Most  practical  skidding  areas  will  have  no  more 
than  two  or  three  different  log  densities.  However,  a 
theoretical  background  should  provide  for  the  gen- 
eral case  and  then  be  adapted  to  specific  instances. 
First  consider  the  general  case  of  variable  log  density. 

Assume  a  skidding  area,  A,  is  made  of  n  subareas, 
A,  of  differing  log  density,  (fig.  5).  The  number  of 
turns  required  to  skid  all  logs  in  area  Aj  is  the  number 
of  logs  in  A.  divided  by  the  average  number  of  logs 
carried  per  turn. 


Pt 

where      tj  =  number  of  turns  required  to  skid  all 
logs  in  area,  Aj. 
1.  —  number  of  logs  in  area,  Aj. 
p^.  =  average  or  estimated  number  of  pieces 
(logs)  carried  per  turn,  and  is  assumed 
constant  in  all  areas. 


The  number  of  logs  in  Aj  is  the  density  per  acre  (or 
other  unit  of  area)  muliplied  by  the  area. 

where:  <$•  =  density  in  logs  per  unit  area  (acre). 

Total  distance  traveled  by  skidding  all  logs  in  area 
Aj  is  the  number  of  turns  multiplied  by  the  AYD  of 
the  subarea,  the  distance  between  the  landing  and 
polar  centroid  of  Aj. 

dj  =  tj  AYD|  =  V, 

where  dj=  distance  traveled  skidding  all  logs  in 
area  Aj. 

r.  =  AYDj  =  distance  from  the  landing  to  the  polar 
centroid  of  area,  Aj. 

Total  skidding  distance  for  all  the  subareas,  Aj,  is 
found  by  summing  dj  for  all  i  areas. 

n 

D=  .2dj. 
i=1  ' 
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Similarly,  total  turns  required  to  skid  all  logs  from  the 
entire  area  A  are 

n 

T  -    X  t, 

1=1  ' 

Average  skid  distance  is  total  distance  travelled  skid- 
ding divided  by  total  turns, 


where  S  and  r  are  functions  of  the  area  involved.  To 
summarize,  the  general  formula  for  average  skid 
distance  depends  on  the  areal  variation  of  log  density. 
S,  and  on  the  distance,  r,  from  deck  to  log  location. 
Shape  of  the  area  is  not  a  factor  in  this  generalized 
formula  but  is  considered  when  applying  the  method 
to  a  specific  skidding  example,  shown  below. 


Id. 


ASD  =  £  = 


Substituting  all  previously  derived  quantities  gives 


ASD  = 


n 


2t, 

1=1 


Average  Skid  Distance  Formula  — 
Circular  Sector 

To  illustrate  the  effect  of  variable  log  density, 
apply  the  derivation  results  to  a  circle  sector  used  as 
an  example  by  Suddarth  and  Herrick  (1964).  In 
addition,  partition  the  circle  sector  into  concentric 
zones  of  differing  log  density  (fig.  6).  This  configu- 
ration, for  example,  might  approximate  a  wedge 
shaped  cable  logging  setting  with  different  cutting 
densities.  For  a  small  differential  area,  da,  in  any 
zone 

da  =rdrdO. 


ASD  = 


i=1  Pt 
~~ "  i 

i=1  Pt 


For  all  n  zones,  ASD 


n     f*  k7i-  n  R; 
\  J      J      8  r2  dr  dB 

1   «/0  R.  1  1 


2   I       I  ^rdrdO 

1^0         R.  , 


n  Both  summation  and  integration  are  necessary 

^  §j  Aj  Tj  here  because  log  density,  S,  is  not  a  continuous 

^gp  _  i  =1  function  over  the  range  of  values  for  r.  In  fact,  for  this 

n  example,  8  is  considered  constant  within  each  log 

^  8 1  Aj  density  zone,  just  like  a  skidding  area  with  several 

'=1  cutting  intensities. 


Note  that  average  or  estimate  number  of  logs  carried 
per  turn,  p^.,  cancels  and  is  not  a  factor  in  the 
remaining  derivation.  This  is  allowed  only  when 
estimated  number  of  logs  carried  per  turn  can  be 
assumed  to  be  constant,  as  discussed  earlier  in  the 
report. 

If  we  let  the  number  of  areas  A^  in  A  become 
infinitely  large,  we  get  the  following  general  integral 
form  for  average  skid  distance: 


ASD 


J 8  (a)  r  (a)  da 
A 


8  (a)  da 


zone 
n-i 


y\         dr  \ 

zone  \  rde         dQ  \ 

i  Sil       82  1 

1    $ri-i|  1 

Rn-2    Rn-i  Rn 


Figure  6.  —Circular  sector  with  n  concentric  log  density 
zones  and  with  radius  R. 
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Proceeding  with  the  evaluation  of  this  formula 
gives 


ASD  - 


ASD  = 


n 

X 

i=1 

8i  . 

Ri-i 

de 

n 

Ri 

de 

i=1 

0 

Ri-i 

I     2  5,(Rf-Rf. 


XkTT 
de 


2    ..(Rf-Rf.,)  f 


de 


ASD  =  £ 


|  8i  (R?-R?-i 
2  8,  (Rf-Rf-i) 


This  is  the  formula  for  average  skid  distance  for  a 
circle  or  circle  sector  with  concentric  zones  of  varying 
log  density.  Note  three  points.  In  figure  6,  RQ  =  0,  but 
this  is  not  a  requirement;  RQ  can  be  greater  than  zero. 
Also,  we  do  not  have  to  include  all  bands;  for 


example,  density,  #2,  could  be  zero.  Finally,  when  all 
log  densities  are  equal,  the  summations  for  average 
skid  distance  reduce  to  ASD  =  2/3R,  the  same  result 
for  a  circular  segment  computed  in  research  cited 
earlier. 

Now  apply  the  general  formula  for  a  circular  sector 
skidding  area  to  a  specific  example,  with  three 
concentric  zones  of  different  log  density,  as  shown 
below: 

Zone  log     Zone  inner      Zone  outer 
Zone,  density,  <5j,  distance  from  distance  from 
1       logs/acre   landing,  Rj..,,     landing,  Rj, 


100 

50 
300 


0 

200 
400 


200 
400 
1,000 


Substituting  these  numbers  into  the  average  skid  dis- 
tance formula  for  a  circle  or  circular  segment  gives  an 
average  skid  distance  of  724  feet.  Average  skid 
distance  is  667  feet  when  variable  log  density  is 
neglected. 

As  noted  earlier,  the  general  integral  form  for 
average  skid  distance  can  theoretically  apply  to  areas 
of  any  shape.  Unfortunately,  direct  integral  evalua- 
tion requires  that  log  densities  and  the  equations 
defining  area  boundaries  be  well  defined,  mathemati- 
cally speaking.  Thus,  evaluation  of  definite  integrals 
is  relatively  easy  for  circles,  circular  segments,  and 
rectangles.  However,  integrals  for  irregular  shapes 
are  difficult,  if  not  impossible,  to  evaluate  and  are 
more  efficiently  evaluated  by  approximate  numerical 
methods  as  illustrated  in  this  report. 
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